Magnetic field effects on two-leg Heisenberg antiferromagnetic ladders: 

Thermodynamic properties 
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Using the recently developed transfer-matrix renormalization group method, we have studied the 
thermodynamic properties of two-leg antiferromagnetic ladders in the magnetic field. Based on dif- 
ferent behavior of magnetization, we found disordered spin liquid, Luttinger liquid, spin-polarized 
phases and a classical regime depending on magnetic field and temperature. Our calculations in 
Luttinger liquid regime suggest that both the divergence of the NMR relaxation rate and the anoma- 
lous specific heat behavior observed on Cu2(C5Hi2N2)2Cl4 are due to quasi-one-dimensional efi'ect 
rather than three-dimensional ordering. 



PACS numbers: 75.10.Jm, 75.40. Cx, 75.40.Mg 

Recently spin ladders have been at the focus of inten- 
sive research towards understanding the spin-1/2 Heisen- 
berg antiferromagnets in one and two dimensions llj-^. 
Experimentally, several classes of materials like SrCu2 03, 
La6Ca8Cu2404i and Cu2(C5Hi2N2)2Cl4 (CuHpCl) have 
been found whose properties can be well described by the 
two-leg Heisenberg antiferromagnetic ladder (THAFL) 
model ^^- For inorganic oxides, the spin gap A was 
found around hOOK jO). So only the low-energy part of 
spectrum can be explored by measuring spin susceptibil- 
ity, NMR relaxation and neutron scattering |||,0,|l2l. On 
the other hand, the organo-metallic compound CuHpCl 
exhibits a very small spin gap Awlliir lO] which allows 
a full investigation of the spectrum by applying a mag- 
netic field (MP). Chaboussant et al. have shown that the 
NMR relaxation rate exhibits substantially different be- 
havior for different ranges of MF in the low temperature 
limit. On this basis these authors proposed a magnetic 
phase diagram 11^. Moreover, the specific heat mea- 
surements show anomalous behavior when the spin gap 
is suppressed by the MF [|l2|,^. Theoretically, some of 
the MF effects on THAFL were discussed by using ex- 
act diagonalization, bosonization, conformal field theory 
and non-linear cr-model approaches [|3|-|9[. In this pa- 
per, we perform the first calculation of the phase diagram 
(more precisely crossover lines between different regimes) 
using the newly developed transfer-matrix renormaliza- 
tion group (TMRG) technique [g^l • Our findings suggest 
the observed divergence of the NMR rate [|ll| and the 
anomalous specific heat behavior ||lj,|l^ are due to the 
MF effects on THAFL i.e., quasi-lD effects rather than 
3D field-induced ordering when H > A. 

The Hamiltonian for the THAFL in our studies reads: 
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plings, respectively. To confront the experimental find- 
ings on CuHpCl, we set J||= 1, J±/J^^= 5.28. 

The TMRG technique we adopt here is implemented 
in the thermodynamic limit and can be used to evaluate 
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very accurately the thermodynamic quantities [2 
well as imaginary time auto-correlation functions 1 22 
at very low-T for quasi-lD systems. Technical aspects 
of this method can be found in Ref. |£l|. In our cal- 
culations, the number of kept optimal states m = 200, 
while the width of the imaginary time slice e = 0.05 are 
used in most cases. We have also used different m and 
e to verify the accuracy of calculations. The physical 
quantities presented below are usually calculated down 
to T < 0.02 (in units of J\\). The lowest temperature 
reached is T = 0.005. The relative errors, being different 
for different quantities, are usually much less than, at 
most about, one percent for derivative quantities at very 
low temperatures. 

We first determine the spin gap by fitting the spin sus- 
ceptibility Xi using the asymptotic formula proposed in 
|||: x=^e-'^/'^/%/T, T ^ 0, based on the quadratic 
dispersion with a gap A for the single magnon branch. 
Fitting numerical results x in the range T G [0.168,1] 
p5| , we obtain A=4.385 which is very close to the value 
A=4.382 obtained using the T=0 DMRG method (with 
250 states kept and extrapolated to infinite size). 
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where S„^i denotes a S= 1/2 spin operator at the i-th site 
of the n-th chain. Jim are the intra- and inter-chain cou- 




FIG. 1. Magnetization versus H for different T 



In Fig. |l|, we show the magnetization curves for dif- 
ferent values of T € [0.02, 10]. For quasi- ID spin-gapped 
systems, it is well-known that at T = 0, M = for H< 
Hc,=A; M = 1 (in units of S) for H>Hc^=Jj_ + 2J|| 
@ and < M < 1 for iJ e {Hc,,Hc^). When T^O, 
M is nonzero for any H. However, the critical behav- 
ior of the magnetization in the vicinity of Hci and iJc'2 
can only be seen at very low temperatures. The be- 
havior of M{H) is elucidated in Fig. 0. First consider 
the T=0 case in (a). We have calculated M(r=0) at 
H = 7.275, 7.25, 7.125, 700 and 4.4, 4.5, 4.625, 4.75 for the 
upper and lower critical points, respectively. The calcula- 
tions for M{T) were done with m=256 down to rR:!0.005 
(needed) for extrapolation to T=0 limit. Then fitting 
M(T=0) at these H gives the following asymptotic form: 



r 0.380 VS-^^ for i/ = i/+ 
^^^^^"\ 1-0.431 /Hc~l? iorH = H^^, ^^> 

in agreement with universal square-root singularities of 
magnetization in gapped systems |£6[ (see also ]15| , p7| ). 
Independently, we obtain Hci = 4.3823 which is even 
more accurate than the value obtained from fitting xC^)- 
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FIG. 2. (a) Solid lines for the asymptotic behavior Eq. (|2|). 
Symbols x(+) for T^O at H^ H^^ (^Ca) '^^'^ up(down) tri- 
angles for T = 0.02 showing the deviation from Eq. (m). (b) 
dM/dH vs H (in legends) for T<To; (c) For 7o< T< Tc; (d) 
For T>Tc- In (c) and (d), a constant a{T) is subtracted and 
derivatives are amplified as seen in parentheses. 

When T>0, depending on the behavior of dM/dH, 
there are three different cases: 1) dA//diJ has a two-peak 
structure shown in Fig. |(b) for r<ro=0.59+*'°'' (posi- 
tivc(negative) numbers in super (sub)-scripts are bounds 
of errors pa]), similar to the r=0 case. 2) It has a 
single peak structure at H=/=0 in Fig. 0(c) for Tq < 
r<rc=4.86_o.o5- 3) It reaches a maximum exactly at 
ff = for T>Tc in Fig. |(d). Suppose 7 is the coeffi- 
cient of the cubic term in the low-H expansion of mag- 
netization, Tc is given by ^{Tc) = 0. 

Based on the above observations, we can construct a 



magnetic phase diagram as shown in Fig. |[ Strictly 
speaking, quantum phase transitions take place only 
at i?Ci,2 foi' T — 0, and "phase boundaries" are just 
crossover lines for T > 0. At T = 0: 1) as i? < i^Ci = ^, 
the band edge of the continuum has S^^^^i — 1, with an 
effective gap AS*- = Hd — H. The ground state is a 
disordered spin liquid and thermodynamic quantities de- 
cay exponentially at low T; 2) As Hc^ < H < Hc2 , the 
gap vanishes and we find a range of linear in T depen- 
dence for the specific heat and finite values for the sus- 
ceptibility, which is characteristic for the Luttinger liquid 
(LL); 3) The ground state becomes fully polarized when 
H > Hc2 ■ There thermodynamic quantities again decay 
exponentially with A^ff = H — Hc^ at low T. When 
T > 0; the LL regime shrinks gradually and disappears 
at T = To and H = Hm, beyond which the system "for- 
gets" about Hci , Hc2 and A ]g3 ■ The other two phases 
continue to exist until T = Tc, and one finds a classical 
regime for T > Tc- 
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FIG. 3. Magnetic phase diagram: the dots fitted by lines as 
phase boundary in-between the LL, the spin-polarized phase 
and the disordered spin liquid, indicate values of H and T 
maximizing dA4/dH. 

Compared with the phase diagram proposed on the 
basis of 1/Ti measurements ]Q, the major difference is 
the absence of the "quantum critical" phase in Fig. ||. 
We note that the requirement for exhibiting the univer- 
sal "quantum critical" behavior J ^ A ]|2^ is not sat- 
isfied in our case. In addition, for the classical regime, 
we found Tc = 1.109_o.oiiA instead of Tc = A. In 
Ref. ] pd| , the phase boundaries for the quantum critical 
phase with two gapped phases are given as the onset of 
the exponential behavior in l/Ti at T = A^jf . (Conse- 
quently, one obtains H = Hci 2 at T = and Tc = A 
from ASf = at iJ = 0.) The divergence of l/Ti pre- 
sumably disappears at the boundary of the LL regime. 
In fact, 1/Ti has contributions from magnon scattering 
with momentum transfer of both q — and q — n. The 
former process corresponds to a larger gap in the contin- 
uum, but contributes substantially to the relaxation [ p9[ . 
In our calculations the critical behavior of M{H) defin- 
ing the quantum phase transitions is identified directly 



at T = 0. When T > 0, a straightforward extension 
of this definition gives rise to all regimes except for the 
"quantum critical" phase. We should also mention that 
when Jj_ -^ 0, one has A ^ 0,To ^ Tq and Hm — > 0. 

Now elaborate more on the temperature dependence 
of M{T) for various given H as shown in Fig. H. Con- 
sider a cooling process. For T > Tc ^ A, M mono- 
tonically but slowly increases for all H, whereas it can 
change non-monotonically depending upon the values of 
H for T < Tc- When H>Hc2, M continues to increase 
and saturates exponentially with A^.r. However, when 
H<Hci , M first goes up, then down and finally decays to 
zero exponentially with AS*. When Hci < H < Hc^, 
there is always a maximum at T 7^ 0, close to the bound- 
ary between the polarized phase and LL for H > Hm ~ 6 
Ipsi , while separating the former from the disordered spin 
liquid phase elsewhere. There is also a minimum for LL. 
The positions of minima are at T = for H > H„i , oth- 
erwise they are close to the boundary of the LL regime. 
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FIG. 4. Magnetization vs. temperature: curves for H = 9, 
8, 7.415, 7.28, 7.125, 7, 6.75, 6.25, 6, 5.75, 5.5, 5, 4.625, 4.5, 
4.385, 4.125, 3.5 from top to bottom. Symbols A (v) denote 
maxima (minima) at each H. Extra points added between 
curves are for maxima at 7.25, 5.875, 5.812, 2, 1, 0.005 and 
minima at 5.875, 5.812. 

The origin of this nontrivial temperature dependence, 
in particular, the presence of minima and maxima at 
low temperatures, is not fully understood. One possi- 
ble interpretation is due to the excitation spectrum in 
the presence of MF. S^^^^i is a good quantum number for 
the Hamiltonian (^ ) and thus different energy bands are 
shifted by —S^^^^iH. When H < Hc2, the ground state is 
not a fully polarized state(FPS), but the low- lying exci- 
tations correspond to positive S^^^^^, and the maximum, 
roughly speaking, corresponds to the "maximally polar- 
ized" state. This is true so far H>Hm- On the other 
hand, the minimum at T^Q originates from low-lying 
states in Si^i=Q subspace. Those states intersect with 
FPS at H — H„i- Notably, Hm also corresponds to an 
extrapolation of the boundary between the two gapped 
phases in the phase diagram ( Fig.^). It is also curious 
to note that in Fig. m the curves are roughly symmetric 



w.r.t. M = 0.5, if we focus on the low temperature part. 
This reflects the particle-hole symmetry of the problem 
in the fermion representation ||l^,|ll[ . 

We now turn to the specific heat which, similar to 
M{T), shows different behavior depending on H. When 
H < Hci in Fig. ^(a), C^ has a single peak structure as 
expected. The MF reduces ASr, and changes dramati- 
cally the line shape near T = 0.5, as i/ ^ H^ . This is a 
signature of approaching the quantum critical point p8[ . 
When i?Ci< H< Hc^ in Fig. |(b) and (d), a second 
peak at low T is developed exhibiting the LL behavior. 
Linear-T dependence is shown in the insets of (b) and 
(d). Moreover, aX H = H^ , the cusp still remains and 
at H ~ H^ a shoulder emerges. When H > Hc2 in 
Fig. |g(c), the shoulder can still be seen for H — H^ , al- 
though the second peak vanishes. At low-T, Cy decreases 
exponentially with A^, ,. We note the cusps and shoul- 
ders appear outside the LL regime but at the vicinity of 
its boundary. For those H at which a larger second peak 
shows up, the local minima are also located outside the 
LL regime, but nearby. 




FIG. 5. Specific heat at various H (legends): For H < Hci 
in (a); Hci < H < Hc^ in (b) and (d); H > Hc^ in (c). Insets 
for low-T behavior: Cv in (a)-(d) and Cv/T in (b) and (d). 
Triangles denote maximum Cv 

It is also instructive to see the MF effects on the maxi- 
mum specific heat C™"^ and corresponding temperature 
Tmax- As seen in Fig. |^, when H is apphed, C™*^^ 
first declines as a gradual response to the splitting. At 
H = Hm, because of the crossing of two kinds of en- 
ergy states discussed above, C™"^ arrives at a minimum. 
Moreover, we surprisingly found that the curvature of 
Tmax changes its sign at Hci , while it reaches a maxi- 
mum at H = Hc2- These interesting features show an 
intrinsic aspect of the MF effects as demonstrated con- 
sistently by the magnetization and the specific heat. 

Finally we discuss the bearings of our numerical re- 
sults on experimental findings of a diverging NMR re- 
laxation rate and peculiar specific heat behavior. Re- 



cently, Chaboussant et al. found that the NMR rate 
anomalously increases for Hci<H<Hc2 when T de- 
creases 1 11 . These authors attributed this anomalous 
increase to quasi-lD behavior. Our results show that 
it is indeed an intrinsic MF effect on the spin ladders. 
In fact, the increasing of the NMR rate starts already 
at the vicinity of the LL regime which is bounded by 
T^l.GK. This should share the same physical origin as 
the occurrence of the cusps and shoulders for C^ at the 
vicinity of the LL regime. On the other hand, Hammar 
et al. 1 12 1 have measured Cy up to H = 9T, which is 
about (5.10 ± 0.17) J||, taking into account the difference 
between the experimental and numerical results for the 
gap due to other interactions [ |l4|Jll| . As seen in Fig. 3 of 
Ref. Pi 
in Fig 



the overall feature is consistent with our results 
|(a) and (b), e.g. the shift of C™*^"^ and T,nax as 
well as the abrupt change for H=6.6T. The development 
of the shoulder and the second peak at low temperatures 
is clearly seen in more recent measurements |13| in full 
agreement with our calculations. When Hci<H<Hc2, 
we notice that a narrow subpeak emerges from the second 
peak at lower temperature (see Fig. 2(b) of Ref. p3[). 
The subpeak might indicate the on-set of the 3D effects 
[ fi2| , P^ , while the second peak represents the magnetic ef- 
fects on the THAFL. On contrary, the NMR rate changes 
smoothly versus temperature in the LL regime. There- 
fore, the anomalous behavior of the NMR rate results 
from the magnetic effects of THAFL as a characteristic 
feature of quasi-lD gap systems. 



0.9 


\ 

\ 

\ 


/ 
/ 


■ k 

\ 


-^_ 


0.7 


\ 


/ 

/ 
/ 




\ 






0.5 


J 


V 


A-^^ 
tf 


/ - 



H„,5 H„ H„^ 10 

H 

FIG. 6. The maximum specific heat C™"^ (to the left) and 
the corresponding temperature Tmax (to the right) versus H. 

In conclusion, we have proposed a magnetic phase di- 
agram for the two-leg ladders. We emphasize that most 
of the striking MF effects show up in the LL regime, giv- 
ing rise to the divergence of the NMR rate and anoma- 
lous specific heat observed in experiments. Moreover, 
this magnetic phase diagram is generically valid for other 
spin-gapped systems as well and the results on M(T) 
should also shed some light on other quasi-one dimen- 
sional fermion systems which involve either a charge gap 
or a band gap, with H replaced by the chemical potential. 
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G. Chaboussant for helpful correspondence. 
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